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Abstract We present results for the optical absorp- 
tion spectra of small-diameter single-wall carbon and 
boron nitride nanotubes obtained by ab initio calcula- 
tions in the framework of time-dependent density func- 
tional theory. We compare the results with those ob- 
tained for the corresponding layered structures, i.e. the 
graphene and hexagonal BN sheets. In particular, we 
focus on the role of depolarization effects, anisotropics 
and interactions in the excited states. We show that al- 
ready the random phase approximation reproduces well 
the main features of the spectra when crystal local field 
effects are correctly included, and discuss to which extent 
the calculations can be further simplified by extrapolat- 
ing results obtained for the layered systems to results ex- 
pected for the tubes. The present results are relevant for 
the interpretation of data obtained by recent experimen- 
tal tools for nanotube characterization such as optical 
and fluorescence spectroscopies as well as polarized res- 
onant Raman scattering spectroscopy. We also address 
electron energy loss spectra in the small-q momentum 
transfer limit. In this case, the interlayer and intertube 
interactions play an enhanced role with respect to opti- 
cal spectroscopy. 
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1 Introduction 

The science of nanostructures is one of the fields of grow- 
ing interest in materials science. Nanotubes, in partic- 
ular, are of both fundamental and technological impor- 
tance; being quasi-one-dimensional (ID) structures, they 
possess a number of exceptional properties. While the 
peculiar electronic structure - metallic versus semicon- 
ducting behavior ~ of carbon nanotubes depends sen- 
sitively on the diameter and the chirality, i.e., on the 
way the graphene sheet is wrapped up into a tube 



boron nitride (BN) tubes display a more uniform behav- 
ior with a wide band-gap (larger than 4 eV) , almost in- 
dependent of diameter and chirality |2E] ■ The potential 
applications of nanotubes in nanodevices are numerous 
U : super-tough nanotube fibres |S] , gas sensors [S] , field 
effect displays [7], and electromechanical devices [S], to 
name only a few. One of the most spectacular exam- 
ples is the realization of field effect transistors both with 
carbon nanotubes and with BN- nanotubes |llj . 

Since carbon and BN nanotubes are routinely pro- 
duced in gram quantities, the challenge now consists in 
having fast experimental tools for characterization of 
nanotube samples and, if possible, single isolated nan- 
otubes. Optical spectroscopic techniques provide us with 
this tool jl2lli;illl4lll5lll()lll7| . The final goal is to find a 
unique mapping of the measured electronic and vibra- 
tional properties onto the tube indices (n,m). To that 
end, the electronic structure and dielectric properties of 
the tubes are two key areas to study. One possible spec- 
troscopic method is optical absorption spectroscopy with 
direct excitations from occupied to unoccupied states. 
For carbon tubes, the energy difference En between cor- 
responding occupied and unoccupied van Hove singular- 
ities (VHSs) in the 1-dimensional electronic density of 
states (DOS) is approximately inversely proportional to 
the tube diameter d. In resonant Raman spectroscopy I14j 
and scanning tunneling spectroscopy (13| , this scaling is 
employed for the determination of tube diameters (and 
(n,m) indexes in resonant Raman jl8|'l. A recent ex- 
ample is the fiuorescence spectroscopy of single carbon 
tubes in aqueous solution, where E22 is probed through 
the frequency of the excitation laser and En is probed 
simultaneously through the frequency of the emitted flu- 
orescent light ^5]. Also optical absorption spectroscopy 
of nanotubes in aqueous solution [W allows for the spec- 
tral resolution of peaks that can be associated with dis- 
tinct types of nanotubes. For the distance between the 
first VHSs in semiconducting carbon tubes, a simple 
TT-electron tight-binding fit yields the relation En = 



2 



A.G. Marinopoulos, L. Wirtz, et al 



£'22/2 — 2ac^clo/ d, where ac^c is the distance be- 
tween nearest-neighbor carbon atoms. The value for the 
hopping matrix element 70 varies between 2.4 eV and 2.9 
eV, depending on the experimental context in which it is 
used. This fact is a clear indication that the above rela- 
tion gives only qualitative and not quantitative informa- 
tion on the tube diameter and/or chirality. Furthermore, 
for small-diameter tubes the band structure completely 
changes with respect to the band structure of large di- 
ameter tubes, including a reordering of the VHSs in the 
density of states and displaying fine structure beyond 
the first and second VHSs ,19 2Q . This structure sensi- 
tively depends on the tube indices and may be probed 
by optical absorption spectroscopy over a wider energy 
range, possibly extending into the UV regime. 

The spectroscopical characterization of macroscopic 
tube samples is made difficult by the fact that in the 
bulk solid and even in bundles the tubes are not per- 
fectly aligned and do not have a well defined diameter 
and helicity, rendering, in the case of carbon tubes, a 
random mixture of semiconducting and metallic tubes 
PT) . Additionally, the role of intertube interaction in the 
spectra needs to be addressed since the tubes are close 
packed and could interact with each other via long-range 
forces induced by the excitations. Only very recently op- 
tical absorption spectra were reported for aligned single- 
wall carbon tubes of a very narrow diameter distribution 
(4±0.2 A) grown inside the channels of a zeolite matrix 
|22j . Geometric arguments predict three possible tubu- 
lar helicities for the range of diameters around 4 A: the 
armchair (3,3), the zig-zag (5,0) and the chiral (4,2) con- 
figuration. Therefore, this particular case serves as an 
important case study where a direct comparison between 
experimental data and theoretical calculations becomes 
possible. Indeed in Ref. T^" we have shown the relevance 
of a first-principles calculation for these small-diameter 
carbon nanotubes by reproducing the polarization de- 
pendence of the measured optical spectra. 

At present, the dielectric response of tubes in the 
frequency range of the electronic interband transitions 
and the collective excitations (plasmons) of the valence 
electrons (up to 50 eV) is not well understood since 
previous studies focused on the low-energy regime and 
excitations |23II24| . For higher frequencies, up to now, 
the majority of the theoretical studies of the response, 
besides model calculations are mostly limited to 

summing over independent band-to-band transitions ob- 
tained within the semi-empirical tight-binding method 
or the density-functional theory (DFT) framework 
|22II27","28" , or to calculations of the joint density of states 
(JDOS) 29 which use the bandstructure with no ex- 
plicit evaluation of the transition matrix elements. As 
we will show below, these approximations are not suf- 
ficient for a full interpretation of the experiments. This 
shortcoming is not due to the quality of the bandstruc- 
ture calculation itself but instead due to the neglect 
of the induced microscopic components in the response 



to the external field, the local field effects (LFE) [T!?j . 
These effects strongly modify the total response for cer- 
tain polarizations of the external perturbation. Also, in- 
duced exchange and correlation (XC) components ob- 
tained beyond the random-phase approximation (RPA) 
might contribute .23 , 30 . Therefore, important questions 
concerning the electron interaction, excitations and screen- 
ing still remain unanswered. Our approach here is to 
determine the spectra by ab-initio calculations incorpo- 
rating important ingredients of the electron response as 
in previous works [T^ll^U] . 

The present paper is organized as follows: after an 
exposition of the theoretical framework in section |21 we 
present in section|21 ab initio calculations of optical spec- 
tra of small single-wall carbon and boron nitride nan- 
otubes. We identify the infiuence of crystal LFE, XC 
effects and intertube interaction in the spectra. We also 
point out the similarities and differences between carbon 
and BN structures which are related to similarities and 
differences in the respective electronic band structures. 
A certain number of comparisons with experiment allows 
us to verify that the chosen ab initio approach is an im- 
portant improvement and well suited for the description 
of the spectra of these systems. In section 0] we anal- 
yse the spectra of the building blocks of the tubes, i.e. 
the graphene and hexagonal BN sheets (layers) including 
the optical spectrum of graphite. This gives information 
about the the interaction between objects (sheet-sheet) 
in the excited state, even when these objects can be con- 
sidered to be isolated in the ground state. We also dis- 
cuss to which extent the dielectric response of the tubes 
can be understood in terms of the response of the sheets. 
This will help to understand whether some aspects of the 
response of the tubes are inherent to the sheets, and if 
so, they could also be observed in other systems of more 
practical interest, e.g. samples comprising a mixture of 
tubes with different diameters and orientation or multi- 
wall tubes. In sectional we present results for electron 
energy loss in order to have some additional validation 
from existing experimental data, and because a compari- 
son between optical and loss spectra allows a deeper dis- 
cussion and understanding of interaction effects. Finally, 
we conclude with sectional with an overall discussion of 
the results. 



2 Theoretical Framework 

Ab initio calculations in the framework of Density Func- 
tional Theory (DFT) have yielded high-quality results 
for a large variety of systems: from molecules to pe- 
riodic solids and structural defects j3I]. These results 
are however mostly limited to quantities related to the 
electronic ground state, whereas additional phenomena 
that occur in the excited state are not correctly de- 
scribed "37. In particular, the self-consistency between 
the total perturbing potential and the charge response 
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induces Hartree and exchange-correlation potentials that 
have to be dealt with. The former give rise to the so- 
called local field effects, whereas the latter can lead, for 
example, to excitonic effects. Today, in the solid state ab 
initio framework two main approaches can include both 
LFE and XC effects piQ; and can be therefore suitable 
for the description of electronic excitations in nanostruc- 
tures. First, Green's functions approaches within many- 
body perturbation theory: here one adds self-energy cor- 
rections to the DPT Kohn-Sham bandstructure and the 
electron-hole interaction is included via the solution of 
the Bethe-Salpeter equation 30 . This approach has given 
excellent results for various bulk and finite systems, but 
it is computationally very cumbersome and not ready yet 
to be applied systematically to more complex systems. 
The time-dependent Hartree-Fock (TDHF) method rep- 
resents a certain level of approximation within this ap- 
proach, where correlation (i.e. screening) in the self en- 
ergy is neglected. LFE, on the other hand, are still re- 
tained. TDHF has been employed, up to now in a semi- 
empirical way, to calculate optical spectra, including ex- 
citons, of carbon nanotubes [221 and also for the 4 A- 
diameter ones |33| . However, the obtained absorption 
peak assignments in the latter case were in disagreement 
with predictions based on the operative dipole selection 
rules for the specific tubular space groups j22ll27| . This 
suggests that in such systems with important metallic 
character a neglect of screening may lead to problems. 

The second way to calculate spectra is given by the 
time-dependent DFT (TDDFT) gl], where all many- 
body effects are embodied in the exchange-correlation 
potential and kernel. It is in practice an approximate way 
but always treats the variations of the Hartree poten- 
tial exactly (right as in both the Green's functions and 
TDHF approaches). This approach, using the (static) 
adiabatic local density approximation (TDLDA) for the 
XC effects, or even completely neglecting the density 
variations of the XC potential (RPA, on the basis of 
an LDA bandstructure) has been applied successfully to 
many finite and infinite systems j3U| . In particular, ex- 
cellent results for the loss spectra of graphite have been 
obtained using this approximation (35) . This is the ap- 
proach we follow in the present work. 

As a first step in our TDDFT approach, we deter- 
mined the electronic ground state of the systems. The 
Kohn-Sham single-particle equations were solved self- 
consistently in the LDA for exchange and correlation ^SB" . 
For the description of the valence-core interaction we 
have used norm-conserving pseudopotentials which were 
generated from free-atom all-electron calculations |37) . 
For the "isolated" geometries of the tubes and of the 
sheets we calculated periodic arrays with a large distance 
between the building blocks in a supercell geometry, in 
order to minimize the mutual interaction. The crystalline 
valence-electron wavefunctions were expanded using a 
plane-wave basis set (with an energy cutoff of 62 Ryd- 
berg for carbon systems, and 50 Rydberg for BN sys- 



tems). A part of the calculations was carried out using 
the ABINIT code j^i- 

The next step is the linear response calculation of 
the independent particle polarizability |39II40| . It in- 
volves a sum over excitations from occupied bands to 
unoccupied bands: 

(n,k|e-'''+'^'n'».k+q)(m,k+q|e'''^+^''nn,k) _ / ^ J, n 

where (m ^ n) means that the indices m and n of the 
first term are exchanged. The result is checked for con- 
vergence with respect to the number of bands |41| and 
the discrete sampling of k-points within the first Bril- 
louin zone. 

Within TDLDA, the full polarizability x is connected 
to x° via|32I 

X = X°+X°(^c + /xc)x, (2) 

where Vc is the bare Coulomb interaction and fxc, the 
so-called exchange-correlation kernel, is the functional 
derivative of the LDA exchange-correlation potential with 
respect to the electron density. By setting fxc to zero, 
exchange and correlation effects in the electron response 
are neglected and one obtains the Random-Phase Ap- 
proximation (RPA). We have carried out our calcula- 
tions in the RPA and also in the TDLDA for certain 
cases. The inverse dielectric function for a periodic sys- 
tem and momentum transfer q is obtained from: 

^g!g' (q) = ^G,G' + ^c(q + G)xG,G'(q) (3) 

with q in the first Brillouin zone and G, G' are reciprocal 
lattice vectors. The absorption spectrum is obtained as 
the imaginary part of the macroscopic dielectric function 

eM(w) = i/£GlG'=o(q ^ 0; ^) (4) 

whereas the loss function for a transferred momentum 
(q-l-G) is given by — Im[£Q^Q(q)]. In inhomogeneous sys- 
tems, e.g. periodic solids, clusters and structural imper- 
fections the inhomogeneity in the electron response gives 
rise to local field effects (LFE) and the Eg.g' cannot 
be considered as being purely a diagonal matrix. There- 
fore, its off-diagonal elements have to be included in the 
matrix inversion. Making the approximation £Q^Q(q) ~ 
l/eG,G(q) corresponds to neglecting the inhomogeneity 
of the response, i.e. to neglecting the LFE. When these 
effects are neglected and, moreover, all transition matrix 
elements in are supposed to be constant, one arrives 
at the widely used approximation that the absorption 
spectrum is proportional to the JDOS, i.e. proportional 
to the sum over interband transitions from occupied (f ) 
to empty (c) states over the Brillouin zone points k, 
E^,c,k'5(eck - et,k - w). 
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3 Optical Absorption Spectrum for carbon and 
BN tubes 

Let us first look at the optical absorption spectra of 
small-diameter carbon nanotubes (all tliree possible he- 
licities), for which experimental results have recently 
become available. The calculations of the spectra jl3| 
were done for the tubes arranged in a hexagonal lat- 
tice with an intertube distance (distance between tube 
walls), equal to Dt = 5.5 A which leads to nearly iso- 
lated tubes. Additionally, for the (3,3) armchair ones we 
repeated the calculations for a solid with a smaller inter- 
tube distance, Dt = 3.2 A, which is close to the interlayer 
distance in graphite (~ 3.4 A). 

The optical absorption spectra for the small diameter 
tubes are diplayed in Fig. ^ (as well as JDOS curves for 
two cases). In the upper panel of the Figure, the JDOS 
divided by the square of the excitation energy, for the 
(3,3) nearly isolated (thin line) and interacting (in the 
solid) (thick line) tubes are shown. It can be seen that af- 
ter an initial very steep decrease (< 1 eV), in both JDOS 
curves there is a gradual increase starting from 2 eV. In 
the case of isolated tubes a sequence of pronounced peak 
structures up to 5 eV is observed: this can be explained 
from the occurrence of direct interband transitions be- 
tween the van Hove singularities of the density of states 
(DOS). These peaks are smeared out in the JDOS of the 
solid where the tubes are strongly interacting 01] . In the 
lower panels of Fig. ^ the calculated absorption spectra 
for light polarizations perpendicular and parallel to the 
tube axis are displayed (as well as the experimental data 
plotted in the inset). The dashed lines denote results 
in the RPA neglecting LFE, continuous lines including 
LFE. For parallel polarization both LFE and LDA-XC 
effects were found to be negligible and, therefore, in this 
case only the RPA results without LFE are presented. 
The fact that both these effects turned out to be negli- 
gible can explain why the peak positions (A, B and C) 
predicted here for the parallel polarization match very 
well the ones found in two recent DFT-RPA studies |27| . 
which completely neglected LFE and XC effects in the 
response. Our calculated peak positions (A, B and C) in 
Fig- Hare in good agreement (to within 0.2 eV) with the 
experimentally observed peak structures for the parallel 
polarization (inset). 

It is also important to note that the present results 
support the same peak-to-helicity assignment as in pre- 
vious works ^22»27^ . Namely, the three observed peaks - 
A, B and C - are due to the (5,0), (4,2) and (3,3) tubes, 
respectively. Such an assignment is also in accordance 
with the dipole selection rules for these tubular space 
groups E11I2ZI. 

Concerning the effect of intertube interaction in the 
absorption we found that it is rather small for this po- 
larization: tube-tube interaction leads only to a broad- 
ening of the main absorption peak (thick solid curve for 
the (3,3) tubes). This can be explained from an increase 



of the energy range of the possible interband transitions 
brought about by the interaction. 

All the discussion of the results up to now seems 
to indicate that LFE may not be needed. Nonetheless, 
for light polarizations perpendicular to the tube axis - 
due to the presence of depolarization - LFE play an 
important role which cannot be ignored. More specif- 
ically, for this polarization the experimental spectrum 
displays vanishing intensity for frequencies up to 4 eV 
(see inset in Fig. dotted curve). Clearly, neither 
the JDOS calculation nor the calculated spectrum with- 
out LFE can capture this effect. Instead, they both pre- 
dict a series of peaks from 2 to 5 eV. The reason that 
RPA-without-LFE fails here is due to the depolarization 
effect |19II23II45| which is created by the induced polar- 
ization charges. The depolarization is only accounted for 
if LFE are included: as it can be seen in Fig. ^ (second 
panel; continuous curve) LFE suppress the low-energy 
absorption peaks and render the tubes almost transpar- 
ent below 5 eV in agreement with the experiment. It 
should not come, therefore, as a surprise that the recent 
DFT-RPA studies |2I| (which did not consider LFE) did 
not reproduce this transparency for the perpendicular 
polarization. 

The TDLDA result, displayed as the dotted curve in 
the second panel, turned out to be qualitatively similar 
to the RPA-with-LFE result: again the low-energy ab- 
sorption peaks are suppressed. This shows that the main 
effect comes from fluctuations of the Hartree-, and not 
from those of the XC-potential (as also obtained for the 
case of graphite |2S|)- LDA-XC effects only cause a small 
(0.3 eV) redshift of the remaining absorption peak to 5.5 
eV, with respect to RPA-with-LFE. This is a character- 
istic behavior of finite systems jSO]. Still there is an open 
question about the role of electron-hole interaction in the 
case of the small band-gap (4,2) tube however the 
good agreement with experiment indicates that its con- 
tribution should be small for the parallel polarization. 

LFE also have a similar drastic impact for the system 
of isolated zig-zag (5,0) tubes (not shown). The manner 
according to which LFE operate may be understood as 
follows: for perpendicular polarization the tubes form an 
assembly of almost isolated, but highly polarizable, ob- 
jects. An applied external field induces hence a local, i.e. 
microscopic, response - the LFE - which strongly weak- 
ens the total perturbation (i.e. it is a "depolarization"). 
The macroscopic response to this weak perturbation is 
only very moderate, because the electrons are localized 
on the tubes. This is totally different from the screening 
in a bulk metal or small-gap semiconductor, where even 
a very weak total perturbation still leads to a strong 
response at low frequencies. For polarization parallel to 
the tube axis, the situation resembles rather this latter 
case, something which explains the absence of LFE for 
this polarization. On the other hand - for the perpendic- 
ular polarization - when the tubes arc interacting in the 
solid (third panel) the depolarization is much weaker: 
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experimentally in a qualitative study of the absorption 
spectra for perpendicular light polarizations. 



2 3 4 5 6 
energy (eV) 

Fig. 1 Calculated optical spectra for the (3,3), (5,0) and 
(4,2) nearly isolated and solid (only for (3,3) case; thick lines) 
of carbon tubes with (continuous lines) and without (dashed 
lines) LFE in the RPA. For comparison the TDLDA result 
for the isolated (3,3) tubes is given by the dotted line in the 
second panel as well as the JDOS curves. The results are pre- 
sented for both perpendicular and parallel to the tube-axis 
polarization. For the (4,2) tube only the result for the par- 
allel polarization without LFE is given, which should be a 
sufficiently good approximation for this polarization, accord- 
ing to our results for the (3,3) and (5,0). The experimental 
data |22| are displayed in the inset. 



the tubes start to absorb (they are no longer transpar- 
ent) because the electronic states start to delocalize and 
the system is now more similar to a bulk metal. This 
different behavior of the response depending on the in- 
tertube distance leads to a very important consequence: 
it suggests that the inter-tube interaction can be detected 
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Fig. 2 Calculated optical absorption spectra for quasi- 
isolated BN nanotubes 47 with (solid line) and without 
(dashed line) LFE in the RPA approximation: a) Compari- 
son of the spatially averaged spectrum of the (3,3) tube with 
the joint density of states divided by the square of the ex- 
citation energy (dotted line); b) spectrum with polarization 
perpendicular to the tube axis for the (3,3) tube; c) spectra 
with polarization parallel to the tube axis for the (3,3), (5,5) 
and (6,0) tubes. 



In Fig. 121 we present optical absorption spectra for 
three different BN nanotubes. First, we compare for the 
case of the BN(3,3) tube the spatially averaged spec- 
trum with the joint density of states (JDOS), divided 
by the square of the transition energy. If local field eff- 
fects are neglected, most peaks of the JDOS are visible 
in the averaged absorption spectra while some peaks are 
suppressed due to small or vanishing oscillator strength 
in Eq. 0J. This demonstrates that the fine structure 
in the spectra is not an artifact of low k-point sam- 
pling but is due to the presence of van-Hove singular- 
ities in the 1-dimensional density of states of the tubes. 
Proper inclusion of LFE leads to a smoothing of the spec- 
trum. However, some fine structure survives and may 
be discernible in high-resolution optical absorption ex- 
periments. The spectrum with perpendicular polariza- 
tion demonstrates that for BN nanotubes, depolariza- 
tion effects play a similarly important role as for carbon 
tubes: Neglecting LFE, the onset of absorption would be 
around 4.5 eV. LFE, however, lead to a redistribution of 
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the oscillator strength to higher energies and render the 
tubes almost transparent up to 8 eV. The band gap of 
BN nanotubes in DFT-LDA is 4 eV and is only weakly 
dependent on radius and chirality |2||S1 • Accordingly, the 
absorption spectra for polarization parallel to the tube 
axis are very similar for the (3,3), (5,5), and (6,0) tubes 
which all display a strong absorption peak around 5 eV 
and a second high peak around 14 eV. We expect these 
structures to be stable also for tubes with larger diam- 
eters. Only the fine structure of the spectra (e.g., the 
absorption peak at 9 eV for the (3,3) tube) depends on 
the details of the band-structure and varies for the dif- 
ferent diameters and chiralities [201 • For the (6,0) tube, 
the first high absorption peak is split and the dominant 
peak shifted towards lower energy. This is a curvature ef- 
fect which leads to a reduction of the band-gap for small 
zigzag BN tubes. 

In the next section we will see to which extent the 
above findings for carbon and BN nanotubes can be un- 
derstood by an analysis based on results for graphene 
and BN sheets and graphite. 



4 Optical Absorption Spectrum for graphite and 
the graphene and BN sheets 

Hexagonal graphite has the ABA Bernal stacking se- 
quence of the graphene sheets. In the present work we 
assumed the experimental lattice parameter ahex and 
{c/a,)hex ratio (2.46 A and 2.73, respectively 0H1)- 

The calculated RPA optical absorption spectra |^ . 
with and without LFE, are shown in Fig. |3 for in-plane 
light polarization (E _L c) and in Fig. 0] (a) for polariza- 
tion parallel to the c axis (E || c). The in-plane spectrum 
is dominated by a very intense peak structure at low fre- 
quencies (up to 5 eV) and also another peak structure 
of broader frequency range which sets in beyond 10 eV 
and has a pronounced peak at 14 eV. The origin of these 
peak structures is due to tt — > tt* and <7 — > ct* interband 
transitions, respectively, according to the earlier inter- 
pretations by Bassani and Paravicini |50| who assumed 
a two-dimensional approximation - no interaction be- 
tween the graphene sheets - and the operative dipole 
selection rules for this polarization. Our calculations of 
the oscillator strength for specific transitions between 
bands in the Brillouin zone (BZ) are in agreement with 
their interpretation. 

LFE arc found to be nearly negligible for this polar- 
ization. This is not surprising since for in-plane polar- 
izations graphite is homogeneous in the long-wavelength 
limit {q — > 0). The general aspects of the spectrum - 
peak positions, their intensity and lineshape - are in 
close agreement with the existing experimental results 
|5T| and the previous all-electron calculation of Ahuja et 
al. [S3] who neglected LFE. 

The absorption spectrum of graphite for the light po- 
larization parallel to the c axis is shown in Fig. 0] (a) . 
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Fig. 3 Absorption spectrum of graphite for E _L c. The star 
symbol denotes the unoccupied electron states. 



It is characterized by a weak intensity in the low fre- 
quency range (0-5 eV) and important peaks in the fre- 
quency range beyond 10 eV. For this polarization the 
bandstructure does not play the exclusive direct role in 
defining the absorption spectrum. Now LFE are very 
important. When LFE are not considered the peak po- 
sitions for this polarization are at 11 and 14 eV as in the 
earlier DFT-RPA calculation However, when LFE 
are included transitions are mixed and the absorption 
spectrum is appreciably modified. The main effect of lo- 
cal fields is to shift the oscillator strength at 10-15 eV to 
higher frequencies. They decrease the intensity of the 11 
eV peak and are responsible for the appearance of the 
16 eV peak in the spectrum. The latter peak is seen in 
experiments as a shoulder suggesting that the inclusion 
of LFE is necessary. 

The non-zero oscillator strength found below 5 eV is 
attributed to the inter-layer interaction which is present 
in the solid. It is also observed experimentally. The dipole 
selection rules for an isolated graphene sheet (i.e. 
complete two-dimensionality) lead to vanishing matrix 
elements and oscillator strength at this frequency range. 
LFE do not have any influence on the lower part of the 
spectrum (less than 10 eV). 

The existing experimental evidence is not conclusive 
for the dielectric function in the 11 eV frequency region. 
The frequency-dependent Im[£M] obtained from electron 
energy loss data [SHl55ll56| displays a very sharp and in- 
tense peak (Im[£M]=10) at 11 eV. On the other hand, on 
the basis of optical measurements the observed max- 
imum at this frequency is of considerably smaller inten- 
sity. The earlier interpretations j57) were based on semi- 
empirical tight-binding bandstructure calculations in the 
two-dimensional approximation (i.e. isolated graphene 
layers) and they predicted peak structure between 13.5 
and 16.5 eV. 

It is therefore also important to understand the ef- 
fect of the inter-layer interaction in the optical response 
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and if this interaction is primarily responsible for the 
occurrence of the intense peak at 11 eV. For this pur- 
pose, we progressively increased the inter-layer spacing 
(doublying and tripling the (c/a);iea; ratio). This yields 
stackings of graphene layers in the unit cell with much 
weaker mutual interaction. The absorption spectra for 
these graphene-sheet geometries are displayed in Fig. 0] 
(b,c). The first observation is that the oscillator strength 
vanishes completely in the frequency region below 10 eV 
in complete accordance with the predictions based on 
the dipole-selection rules for an isolated graphene layer. 
Therefore, the double inter-layer spacing leads to non- 
interacting graphene layers as far as the RPA absorp- 
tion spectrum (at this frequency range) is concerned. 
When LFE are neglected, the peak structure in the 10- 
15 eV range stays intact, now with a smaller intensity 
due to the larger volume. With increasing interlayer sep- 
aration, LFE become progressively more important. The 
shift of oscillator strength induced by LFE is so big that 
the absorption peaks at 10-15 eV are almost completely 
suppressed. These findings demonstrate that both inter- 
layer interaction and LFE influence considerably the in- 
tensity of the absorption peak at 11 eV. Qualitatively, 
the depolarization effects found in the tubes can hence 
be explained by the local field effects observed in the 
graphitic response perpendicular to the graphene sheets. 
Likewise, the effect of intertube interaction in the spec- 
tra for polarizations perpendicular to the tube axis is 
also consistent with the increased propensity of LFE to 
shift oscillator strength at higher frequencies when the 
intersheet distance progressively increases (see Fig.0J. 



We conclude this section with a comparative presen- 
tation of the absorption spectra of the graphene and BN 
sheets using the same inter-sheet distance for both cases, 
shown in Fig. [S] Also we make comparisons with the 
spectra of the corresponding tubes. The spectrum of the 
graphene sheet for the in-plane polarization resembles 
closely the in-plane polarization spectrum of graphite 
(Fig. 121) (except for a scaling due to the change in the 
volume), confirming once more, as in the case of the 
tubes for the polarization parallel to the tube axis, that 
the position of the absorption peaks and their lineshape 
is only weakly influenced by the distance and inter-sheet 
interaction in this case. The main difference between the 
graphene and BN-sheet spectra for the in-plane polar- 
ization is the complete absence of any feature below 4 
eV in the BN spectrum. This is clearly related to the 4 
eV LDA-band gap of BN. The out-of-plane polarization 
spectra of the graphene and BN sheets are remarkably 
similar. Most importantly, both display a transparency 
up to 10 eV and both demonstrate a strong depolar- 
ization effect with shift of oscillator strength to higher 
energies PH) . 

We compare now the calculated sheet- absorption spec- 
tra for in-plane polarization with the corresponding tube- 
absorption spectra for polarization parallel to the tube 
axis (Fig.^for carbon and Fig.[21for BN). In the carbon 
case, the strong absorption feature of the sheet between 
and 4 eV maps onto one or several absorption peaks of 
the tubes in this energy range. The details of this map- 
ping are, however, sensitively dependent on the diameter 
and the chirality of the tubes, since dipole selection rules 
play a strong role in these highly symmetric systems. In 
particular, for tubes with larger diameter than the ones 
calculated in this article, the distance En of the first 
van Hove singularities strongly depends on whether the 
tube is metallic or semiconducting. Therefore, calcula- 
tions on graphite or graphene alone will not be sufficient 
to predict absorption spectra of small-diameter carbon 
nanotubes. For BN, in contrast, the first high absorption 
peak at 5.5 eV maps directly onto a corresponding peak 
in the tubes. Also the high energy absorption feature 
around 14 eV is very similar for the BN sheet and the 
BN tubes. Due to the large band gap of BN, the com- 
parison between sheet and tube spectra is much more 
favourable for BN than for carbon tubes. 

On the other hand, in the case of carbon, the ex- 
trapolation of results from graphite and graphene to the 
tubes is much more straightforward for the electron en- 
ergy loss spectra as discussed in the following section. 

5 EELS spectra in the limit of small q: tubes 
versus the layered structures 

Having seen the strong interaction effects that occur in 
optical spectra for a polarization perpendicular to the 
planes or tube axis, it is interesting to make an excur- 
sion to a different type of spectroscopy, namely electron 
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Fig. 5 Calculated absorption spectra for the graphene and hexagonal BN sheets |58|. 



energy loss (EELS) Although also in this case one 
can measure, like in the absorption experiment, a mo- 
mentum transfer close to zero, there is a substantial dif- 
ference in the definition of the response function: in the 
absorption measurement, one detects the response to the 
total macroscopic field, whereas in EELS the response to 
the external field is reported. Therefore, as pointed out 
above, absorption is linked to the imaginary part of the 
macroscopic dielectric function, but EELS is linked to 
the imaginary part of the inverse of the latter. Math- 
ematically this translates into the fact that the crucial 
response function for EELS is governed by equation (O , 
whereas in the case of absorption one can use slightly dif- 
ferent quantity: the macroscopic dielectric function can 
in fact be rewritten as 

EAiiuj) = 1 - limq^oVciq, G = 0)xG=G'=o(q; uj) (5) 

where Vc{G — 0) is the long-range component of the 
Coulomb potential, and it has a 1/q^ divergence for van- 
ishing q and the quantity x obeys an equation similar to 
(jSJ: X = x" + x"{V'c + fxc)X7 but setting the G = com- 
ponent of Vc to zero (Vc) 30 . Therefore, this seemingly 
tiny difference is responsible for the difference between, 
e.g., the position of the main absorption peaks and that 
of the valence plasmons in solids, and one can also ex- 
pect that it will be crucial when interaction effects on 
the spectra are discussed. In particular, EELS spectra 
should show, due to the presence of this long-range term, 
stronger interaction effects than absorption spectra. 

This is in fact the case, as we will illustrate in the 
following for the total n + a plasmon in graphite. This 
plasmon represents the collective excitation mode of all 



the valence electrons in graphite. Before discussing the 
results we stress here that the tubes cannot be consid- 
ered as completely isolated objects in this calculation 
of the loss function with the 5.5 A intertube distance. 
We refer to the tubes in the latter geometry as distant. 
Fig.Elshows the RPA loss function, —Im{l/eM), for the 
(3,3) tube, in the range 15-35 eV and for a vanishing 
momentum transfer q parallel to the tube axis. For this 
orientation LFE are negligible. A strong shift of the vt+ct 
plasmon from 22 to 28 eV due to intertube interactions in 
the solid can be seen |E3 • The magnitude of this shift re- 
veals a strong dependence of the plasmon position upon 
the intertube distance (hence the average valence elec- 
tron density) essentially following the plasmon-frequency 
dependence in the case of the homogeneous electron gas 

m 

This shows that the tubes respond as homogeneous 
and highly polarizable objects for parallel q in the long- 
wavelength limit (q — s- 0). A direct consequence would 
then be that the atomic arrangement, orientation of bonds 
and hclicity may play a secondary role in the response in 
this frequency and q range. Therefore, the result for the 
n + a plasmon (shown in Fig. EJ would be representa- 
tive of either of the three tubes since all of them - being 
of nearly the same diameter ~ possess the same aver- 
age electron density. This is indeed the case as it can be 
seen in the inset of Fig. El where an almost indistinguish- 
able TT + a plasmon was also obtained for the (5,0) tube. 
Hence, in contrast to an optical absorption experiment, 
small-q loss measurements of the tt + a plasmon cannot 
determine tubular helicities for a given tube diameter. 
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The governing factor for the n + a plasmon must 
be traced to the in-layer graphitic response as it can 
be seen in Fig. where the loss function for graphite 
and graphene is also shown (for in-layer q orientation) 
at comparable (to the tubes) average electron densities 
(to within 10 %) (dashed curves in Fig. EJl. This 
shows that the loss function of the tubes for parallel q in 
this frequency range is governed by the average-density- 
dependent part of the in-layer graphitic response. Similar 
plasmon shifts, therefore, can also be expected in other 
carbon systems with graphene-based structural blocks 
e.g. multiwalled tubes. The present results outline the 
significance of the tt + cr plasmon as a key measurable 
spectroscopic quantity which could gauge the intertube 
distances and interactions in real samples of carbon nan- 
otubes. 




loss energy (eV) 



Fig. 6 Calculated RPA loss function for the carbon (3,3) 
distant tubes (continuous thin line) and solid (continuous 
thick line) of tubes, q is vanishing with an orientation parallel 
to the tube axis. The dashed curves denote the results for 
graphite (thick) and graphene (thin) for small in-layer q j^]. 
The loss function for the (5,0) tube is shown in the inset as 
a dotted curve. LFE are negligible for this q orientation. 

In order to gain a more complete understanding of 
the in-layer graphitic response at small q's and how the 
latter is influenced by the inter-layer interaction we de- 
termined both the loss and dielectric function for various 
graphene-like geometries i.e. varying the interlayer spac- 
ing or, equivalently, the (c/a);iea: ratio. In these calcula- 
tions we also looked at the lower-frequency tt plasmon. 
Bearing in mind the discussion in the previous para- 
graph, these calculations could then serve as benchmarks 
for predicting the position of the tt -I- ct plasmon in nan- 
otubes as a function of the intertube distances for low 
q's. 

The loss and dielectric function for small in-plane 
q (0.22 A^^) is shown in Fig. [7| for graphite and the 
graphene geometries with multiple {c/a)hex ratios. It can 
be seen that the peak positions of both plasmons have 
shifted to lower frequencies when the inter-layer sepa- 
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Fig. 7 Loss and dielectric function for small in-plane q (= 
0.22 A~^) for graphite and the graphene-like geometries with 
multiple (c/a)hej; ratios. Only the results without LFE are 
presented since LFE are negligible in this q range. 

ration is increased. This effect is very pronounced for 
the TT + cr plasmon position. These results reaffirm that 
the total (vr + <t) plasmon is extremely sensitive to the 
inter-layer interaction for small in-plane q's. This can 
be explained as follows: with increasing (c/a)?iea;i i-e. 
interlayer spacing, the system becomes an assembly of 
nearly isolated graphene sheets. This leads to a decrease 
of screening (Re[eM] 1; see Fig. [TJ with the effect 
that — /m[l/e7\/] — > /TO[eA/], namely a coincidence of 
the loss and absorption functions. Since for the in-plane 
polarization the positions of the absorption peaks do not 
change with increasing (c/a)/iea; (see Im[£M] in Fig. [71), 
then the loss function undergoes important changes, in 
particular the -KA-a plasmon. The latter is displaced at a 
much faster rate to lower frequencies towards the 14 eV 
peak of the absorption spectrum whereas the tt plasmon 
peak is rather insensitive to (c/a);iea; since it is already 
located very close to the 0-5 eV peak structure of the 
TT — > TT* transitions in lm\eM\- 

At present, existing measurements [6211 (iHj of tlie loss 
spectra of samples of single- wall carbon nanotubes have 
given a TT -f cr plasmon in the frequency range 21-24 eV 
for small momentum transfer q. Our predicted frequency 
of the TT -t- cr plasmon for the case of the distant tubes 
is within this frequency range (FigEJ. However, a direct 
comparison of the present results for the loss function of 
(3,3) and (5,0) tubes with the measured loss data is not 
straightforward. The difficulty stems from two factors, 
tube diameter and alignment, which have a competing 
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effect on the tt + cr plasmon position: a) Bulk samples 
of single-wall tube material possess a mean diameter of 
14 A IHS], which is considerably larger to the range of 

4 A studied in the present work. Assuming a common 
intertube spacing, then larger diameters would give rise 
to a displacement of the vr + cr plasmon towards lower 
frequencies since the diametrically-opposed wall parts of 
the same tube are facing each other at larger distances, 
i.e. a situation resembling a stacking of graphene sheets 
with larger {c/a)hex ratios (see Fig.[7I). For instance, the 

5 {c/a)hex ratio corresponds to an intersheet separation 
of 16.8 A and the corresponding tt -f cr plasmon peak is 
at 16 eV. b) The alignment of the tubes in the samples is 
not perfect; therefore it is to be expected that also out- 
of-plane excitations of graphitic origin will contribute to 
the response. These excitations [SS] should tend to pro- 
duce a more diffuse shape for the loss spectrum, heavily 
dampening the tt -|- cr plasmon and shifting the observed 
peak towards higher frequencies. 

The type of response described just above can be 
clearly seen in Fig. |S1 which shows the loss function of 
the distant (3,3) tubes for q orientation perpendicular 
to the tube axis. For this orientation in- plane as well 
as out-of-plane graphene excitations contribute to the 
tube response. The latter cause the diffuse shape of the 
loss function (see inset of Fig. (SJ. LFE are now very 
important and the peak position of the tt + cr plasmon 
is at 28 eV. 



collective oscillation (tt -f cr plasmon) around 20 eV. The 
exact position of the peaks depends on the radius and 
chirality of the tubes and on the intertube distance. In 
order to compare with experimental EELS-spectra on 
multi-wall BN tubes p5, an extrapolation to larger- 
diameter tubes is needed and the inter-wall interactions 
have to be taken into account. 



n ' r 




energy (eV) 



Fig. 9 Loss spectra for three different BN nanotubes with 
momentum transfer q ^ Q along the direction of the tube 
axis. The spectra are calculated with an inter-tube distance 
of 7.4 A. Since the spectra with and without LFE are almost 
indistinguishable, only the calculation without LFE is shown. 
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Fig. 8 Calculated RPA loss function for the carbon (3,3) 
distant tubes for small q (0.19 A^^) of orientation perpen- 
dicular to the tube axis. In the inset the loss function of 
graphene |35| for 2 {c/a)hex and for a similar q of orienta- 
tion perpendicular to the sheet is displayed. Continuous and 
dashed curves denote results obtained with and without LFE, 
respectively. 



For completeness, we show in Fig. |^ also the calcu- 
lated EELS spectra for small BN nanotubes with mo- 
mentum transfer g — > along the tube axis. As in the 
case of carbon tubes, two main features are clearly pro- 
nounced: the TT-plasmon at 6-7 eV and the high energy 



6 Discussion and Conclusions 

In this closing section we should like first to comment 
on the validity of RPA and TDLDA for the description 
of the optical spectra of the tubes. It is well known that 
RPA and TDLDA often give excellent results for loss 
spectra [3D]. However, in this work we have seen that 
also very good absorption spectra were obtained for the 
carbon tubes, even at the RPA level, despite the fact 
that both RPA and TDLDA are known to fail badly in 
the description of absorption spectra of many bulk ma- 
terials (silicon and argon being two representative cases) 
|3(J) . Our explanation in this regard is the following: for 
the light polarization parallel to the tube axis, the en- 
suing screening is significant and, therefore, XC effects 
are damped. For the perpendicular polarization and for 
the larger inter- wall distance {Dt — 5.5 A) the tubes be- 
have essentially like isolated systems, where strong can- 
cellations are known to occur between self-energy cor- 
rections and the electron-hole interaction, i.e. between 
XC effects (see Ref.|30|'). The experimental precision is 
then not high enough (also due to the almost vanished 
absorption intensity in the relevant frequency range) to 
discern to which extent XC effects should be better de- 
scribed by approximations beyond the TDLDA. In view 
of these considerations, we can be confident regarding 
the quality of the calculations. On the other hand, the 
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results for the strongly interacting tubes (smaller inter- 
tube distance) should be regarded as qualitative since 
the system becomes then more similar to a solid where 
the cancellations may be more incomplete, and since no 
direct comparison to experiment is possible, at present, 
in this case. 

Still, excitonic effects are not included in the calcula- 
tions. They could play a role in both carbon and BN nan- 
otubes, leading to redistribution of oscillator strength 
and/or appearence of new peaks in the bandgap (bound 
excitons). They might be the reason for the anoma- 
lous E11/E22 ratio measured recently for carbon 
tubes ^ni- We expect a stronger deviation of measured 
optical spectra from theoretical ones in the case of BN 
nanotubes. For this large band gap material, two effects 
will most likely play an important role: Quasi-particle 
corrections will widen the band gap as in the case of 
hexagonal BN, where a GW calculation has demonstrated 
a band-gap increase from 4 eV to 5.5 eV Excitonic 
effects, in contrast, will lead to isolated states in the 
band gap or to an overall reduction of the band-gap. 
To which extent quasi-particle corrections and excitonic 
effects cancel each other is presently not clear. Work 
along these lines is in progress (66j . For the role of quasi- 
particle corrections and excitonic effects in carbon tubes, 
we refer the reader to Ref. |46|. 

In conclusion, the present ab initio calculations of 
the optical absorption of small-diameter carbon and BN 
nanotubes give good agreement with the available ex- 
perimental data. The inclusion of local field effects in 
the response to a perturbation with perpendicular po- 
larization is necessary for a proper description of the 
depolarization effect leading to a suppression of the low- 
energy absorption peaks for both types of tubes. This 
suppression can also explain recent findings in near-field 
Raman microscopy of single- wall carbon nanotubes |67|. 
The proper analysis of the polarization dependence of 
the absorption cross section is very important in order to 
describe the surface enhanced Raman scattering exper- 
iments in isolated carbon nanotubes |18II68) . In carbon 
tubes the position of the first absorption peak strongly 
varies with the tube indices while in BN tubes the first 
peak is determined by the band gap of BN and is there- 
fore mostly independent of (n,rn). For the BN tubes 
some of the fine-structure which distinguishes tubes of 
different chirality is only visible in the UV region which 
gives rise to the hope that this energy regime will be 
probed in the future. 

The intertube interaction was also found to be very 
important. For the carbon tubes this interaction is the 
governing factor which determines the position of the 
higher-frequency t: + a plasmon. This plasmon, hence, 
may prove to be a very useful spectroscopic quantity 
probing intertube interactions and distances in real sam- 
ples. Finally, the corresponding results for the layered 
constituents — graphene and BN sheets — revealed that 
some aspects of the tubular dielectric response can be ex- 



plained even at a quantitative level from the in-layer and 
interlayer response of the sheets. 
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